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I. INTRODUCTION 



Recent developments of the string theory suggest that in a regime of Planck length 
curvature, quantum fluctuation is very large so that string coupling becomes large and 
consequently the fundamental string degrees of freedom are not a weakly coupled good ones 
n. Instead, solitonic degrees of freedom like p-brane or D-p-branes M are more important. 
Therefore it is a very interesting question to ask what is the effect of these new degrees of 
freedom to the space time structure especially whether including these degrees of freedom 
resolve the initial singularity, which is a problem in standard general relativity. 

For the investigation of the p-brane cosmology, the usual low energy effective action 
coming from the beta function of the string worldsheet would not a good starting point. 
So there will be a difference from the string cosmology ||. We need to find the low energy 
effective theory that contains gravity and at the same time reveals the effect of these solitonic 
objects. But those are not known. Therefore one can only guess the answer at this moment. 
The minimum amount of the requirement is that it should give a gravity theory therefore 
it must be a generalization of general relativity. The Brans-Dicke (BD) theory 0] is a 
generic deformation of the general relativity allowing variable gravity coupling. In fact low 
energy theory of the fundamental string contains the Brans-Dicke theory with a fine tuned 
deformation parameter (u=-l). Moreover Duff and et al. |J found that the natural metric 
that couples to the p-brane is the Einstein metric multiplied by certain power of dilaton field. 
In terms of this new metric, the action that gives the p-brane solution becomes Brans-Dicke 
action with definite deformation parameter uo depending on p. Using this action, Rama || 
recently argued that the gas of solitonic p-brane || treated as a perfect fluid type matter 
in a Brans-Dicke theory can resolve the initial singularity without any explicit solution. In 
a previous paper 0, we have studied this model and found several analytic solutions for a 
few values of parameter with which the coupled dilaton-graviton system could be decoupled 
by the simple "completing the square" method. In this paper, we give exact cosmological 
solutions for any Brans-Dicke parameter u and for general equation of state and classify the 
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cosmology of the solutions according to the range of parameters involved. 

The rest of this paper is organized as follows. In section II, we set up the action for the 
p-brane cosmology. In section III, we find an analytic solution for the equation of motion 
and constraint equation for the general case. In section IV and V, we study the cosmology 
of the solution and classify them according to their behavior. In section IV, t as a function 
of the dilaton time r is studied and the behavior of the scale factor a with respect to the 
dilaton time r in the asymptotic region is studied in section V. In section VI, using the 
results in section IV and V, we classify the cosmology into several phases and investigate 
the behavior of the scale factor a as a function of the cosmic time t. 

In section VII, we summarize and conclude with some discussions. 

II. CONSTRUCTION OF THE ACTION WITH THE SOLITONIC MATTER 

We consider the bosonic part of the effective string action and analyse the evolution of 
a D dimensional homogeneous isotropic universe with the solitonic matter included. The 
action is given by 

S = J d D xV=ge^ [TZ - ud^dy] + S m , (1) 

where is the dilaton field and S m is the matter part of the action. In string theory the 
BD parameter u is fixed as —1. In the high curvature regime, the string coupling is also big 
and the solitonic p-brane will be copiously produced since they become light and dominate 
the universe in that regime. Duff et al. || have shown that in terms of metric which couples 
minimally to p-brane(p=d-l), the effective action can be written as Brans-Dicke theory with 
the BD parameter uj given by 

(D-l)(d-2)-d 2 
u ~ (jD _ 2 )(d-2)-cr [) 

In 4-dimension, the BD parameter is given by u = — I for the 0-brane {p = 0) and uj = — | 
for the instanton (p = —1). Let's assume that the gas of solitonic p-brane can be considered 
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as perfect fluid in the Brans-Dicke theory with the equation of state p = 7/?, 7 < 1. Therefore 
our starting point is the equation of the BD theory ||, || 

n, v - 9 -^n = ^ + u{d,<j>d v( j> - ^f(d<p) 2 } 

+{-d„d u (t) + d^d v (j> + gflu V 2 <f) - g^(d<f)) 2 }, 

11 - 2uoV 2 (f) + uo(d(j)) 2 = (3) 

where <fi is the dilaton and T> means a covariant derivative. 1Z is the curvature scalar and 
cosmological metric is given as the following form 

ds D 2 = -jfdt 2 + e 2a(t) 5 ij dx i dx j (i,j = l,2,---,D — 1), (4) 

where e a<yt \= a(t)) is the scale factor and J\f is the (constant) lapse function. Now, we 
assume that all variables are the functions of time only. The curvature scalar flO| in D 
dimension is given by 

n = g 00 n 00 +sr i ii ij , 

g O0 K O o = ^^-[a + a 2 ], 

g^Kv = ^±[S + (D - l)a 2 } (5) 

where at means the time derivative of a. 

The energy-momentum tensor of the solitonic matter is given by 

T^u = V9iiv + + p)U tl U v (6) 

where is the fluid velocity. The hydrostatic equilibrium condition of energy-momentum 
conservation is 

p+(D-l)(p + p)a = 0. (7) 
Using p = 7p, we get the solution 

p = p e^ D - 1 ^ 1+ ^ a . (8) 



The parameters 7 and ui expressed in eq.(l) and eq.(8) are free parameters. Our goal is to 
study how the metric variables change their behavior for various values of 7 and uj. 

If we consider only the time dependence, the action can be brought to the following form 

1 



S = Jdt e^ a - 



{ — (D — 2)(D — l)a 2 + 2(D - l)6uf> + u<j) 2 } 



(9) 



where we eliminated p and p by eq.(7) and eq.(8). The variation over the constant lapse 
function, which is only g 00 , gives a constraint equation. When we set the lapse function M 
to be 1 after varying of the action over the lapse function Af, this constraint equation is the 
equation of motion of the goo component in eq.(5). 

III. ANALYTIC SOLUTION 



Now, we introduce a new time variable r by 

dt = e^-^dr. 



Then the action becomes 
S = jdr 



{{D - 1)kY 2 + pX 2 } - VATpoe 



-IX 



where new variables presented in the action are given by 

k=(D-1)(1- 1 ) 2 (u-u k ) j 

v = 2(1 - 7) (a; - u u ), 
4(D-2 



V = 



K 



-2X = (D — 1)(1 — 7) a — 0, 
Y = a + -X, 

K 

D - 2D7 + 2 7 



UJ,, = 



(D-l)(l- 7 ) 2: 
1 



1-7' 
D - 1 

D-2' 



(10) 



(11) 



(12) 
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The constraint equation is written as 

= (D — 1)kY 2 + fiX 2 + p e~ 2X , (13) 
where p is a positive real constant. The equations of motion are written as 

= Y, 

= x - ^e~ 2X . (14) 
fi 

Note that c<j_i in eq.(12) happens to be the value of the instanton. If u is less than the 
kinetic term of the dilaton has a negative energy in Einstein frame. So we will consider the 
case where cu is larger than oj-\. According to the sign of ac, the types of solutions are very 
different. 

When k is negative, an exact solution becomes 

q 

X = In [- cosh(cr)], 

Y = At + B, (15) 

where c, A, B and g = are arbitrary real constants. Using the constraint equation, we 
determine A in terms of other variables 



A=-, with <J = W-- — = — ; 16 

* V /i 2^1 + cgef 

If ac is zero, then we can obtain a solution of the equations of motion, but it does not satisfy 
the constraint equation. If k is positive, the solution is 

q 

X = In [- | sinh(cr) | ], 

c 

Y = -r + B. (17) 
o 



IV. COSMOLOGY OF THE SOLUTION 

Now, we investigate the relation between the cosmic time t and the dilaton time r. Since 
the solutions of the equations of motion have different forms, we study the behavior of t as 
a function of r case by case. 



A. k < case 



In this region, oj < uo K . 7 is always less than I. We find the relation between t and r 
using eq.(10) 



f - to = / T rfr' 



exp ( - P - 1)7 V -(2+ (jP - 1)7t/ )ln{^cosh(crQ} 

K C 



+(D - 1) 7 £ 



(18) 



where (Z) — 1)7-8 is a constant. This constant can be ignored in the limit r — > ±00. Because 
j| is always positive definite, t is a monotonic function of r. The behavior of a(t) as a 
function of t depends crucially on the relation between t and r. When r goes to ±00, t is 
reduced to 



(19) 



where 



2c 



(D- l)74/l + w 



± {« + (£)- 1)7(1 + 7))} 



(20) 



£>- 1 

We define a new concept for our purpose: t is super-monotonic function of r if it is 
monotonic and t runs entire real line when r does so. When t is super-monotonic function 
of r, the universe evolves from infinite past to infinite future. Otherwise the scale factor 
a(t) has a starting (ending) point at a finite cosmic time tj (tf) which corresponds to initial 
(final) singularity. As a mapping, t maps the real line of r to 

(-00,00) if T_ < < T+, 

(-00, t f ) if T_ < and T + < 0, 

(U, 00) if T_ > and T+ > 0, 

(ti,t f ) if T + <0<T_. 



In the limit r — > ±00, the condition T± < is expressed as 
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(D - l) 7 Wl + u;^— I < =f[« + (D - 1) 7 {1 + w(l - 7)}]. 



(21) 



This inequality is divided into two cases according to the sign of 7. In each case we obtain 
the different region of uo satisfying the condition T± < 0. 



1. 7 > case 

Because we have considered only the case uo > u;_i, so the left hand side in eq.(21) is 
positive definite. To satisfy the inequality T_ < 0, the conditions 

k + (D- 1) 7 {1 + - 7)} > and 
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(D - 1)7^1 + uj^—j 1 < (h + (D - lb{ I + ^( L - ,) [)- (22) 

must be satisfied. If they don't, we know that T_ is positive. The first inequality in eq.(22) 
is reduced to the following inequality 

D — (D — 1)7 

- > «- : = (23) 

It is remarkable that the second inequality in eq.(22) is written as 

(uj - uj )(uj - uj k ) > 0, (24) 
where oj k appeared in the definition of k and 

^0 = -^37 (25) 
is the value of u for the 0-brane. See eq.(2). 

Figure 1. 

As shown in Figure 1, uj k < ujq {uj k > oj ) in the region 0<7<yj(7>^). Therefore 
the solution of eq.(24) becomes 

r n 2 

< < w K tor < 7 < — , 
2 

c<j_i < c<j < cj for 7 > — . (26) 



Combining eq.(23) and eq.(26), we find the region of lo satisfying the condition T_ < as 
the following 

1 2 

w_i < u < u K for < 7 < — , 

2 

c<j_i < < cj for 7 > jj- (27) 

The condition T + < is 



(D - 1)J 1 + uj^- j <-[«+(£>- 1) 7 {1 + w(l - 7 )}]. (28) 



For this, two conditions 



k + (D- 1) 7 {1 + w(l - 7)} < and 

2 



^(D - 1)7^1 + <("[« + (£>- 1)7{1 + w(l " 7)}]) 2 , 

must be satisfied at the same time. In eq.(29), the first inequality gives 



(29) 



(30) 



and the second inequality gives eq.(24) again. Therefore using eq.(26) and eq.(30), we find 
the region of uo satisfying T + < 

cj_! < uj < uj k for < 7 < — - . (31) 



2. 7 < case 
In this case, the condition T_ < is written as 



(D - 1) I 7 I Jl + cu^-^ >-[k+(D- 1) 7 {1 + w(l - 7)}]. (32) 



For this, we need 



k + (D- 1) 7 {1 + u(l - 7)} > or 
(D - l)Jl + u Ezl\ >(_[« + (£>_ 1) 7 {1 + c(l - 7 )}]) 2 • (33) 



Eq.(33) can be simplified as 



UJ-oo < UJ OT uj < uj < 0J K . 



(34) 



Thus the solution, which is the sum of two regions in eq.(34), is reduced to 

uj < uj < uj k for 7 < 0. 

This solution includes the region of the first inequality of eq.(34). 
Similarly, the condition T + < is written as 



(D - 1) | 7 | J 1 + uj^-^- >k + (D- 1) 7 {1 + uj(1 - 7)}, 



which gives 



k + (D - 1) 7 {1 + uj(1 - 7)} < or 

2 



{{D- 1)7^1 + >(«+(^-1)7{1 + w(1-7)}) S 



(35) 



(36) 



(37) 



The solution of these can be written as u < U-co or uj < uj < u K . From these, the region 
satisfying T + < is 



uj_i < uj < uj k for 7 < 0. 



(38) 



B. k > case 



Now we consider positive k, which means 



(39) 



Since the solution X{r) has a singularity at r = 0, we have to treat carefully the behavior 
of t near r = 0. The relation between t and r is given by 



t-t 



= / V exp [l^zillV - (2 + l n{ ^ I sinh(cr0 | } 

J T K C 

+(D-1) 7 S . 



(40) 



In the limit r — > 0, the above equation is reduced to 



t-tn 



sign(r) 



q - J?e (D-i) 7 B 



11-7) _ 



T 



I 1-7? 



(41) 



where 77 = 2 + ^ and To and to are real constants. In case of 77 > 1, t has a singularity 
at r — > 0. In the other case, t has no singularity. So we consider two cases rj < 1 and 77 > 1. 



i. 77 < 1 case 

In this case, t has no singularity at r = 0. So we investigate the behavior of t at r — > ±00 
only. 

i) 7 > case 

In the case k > 0, eq.(40) is reduced to 

t-to^(e T±T -e T ^>), (42) 

where 

2c 



T± = 



K 



(D - 1)J 1 + uj^- I =F {/t + (-D — 1)7(1 + w(l - 7))} 



(43) 



The condition T_ < is written as eq.(28) and gives the solution 

uj < uj-00 and u) > ujq for < 7 < — , 

2 

uj < uj^oq and uj > uj k for 7 > — , (44) 

where we use uj > u K . As shown Figure 1, uj q > cj_ 00 for < 7 < ^ and uj k > uj _oo for 
7 > |j. Therefore there is no solution satisfying the condition T_ < 0. Hence T_ is positive. 

Now we investigate the behavior of t at r — > +00. The condition T + < is written like 
eq.(22). Applying the similar method used in the above analysis, the region of uj satisfying 
T + < is summarized as the following 

ujq < uj for < 7 < — , 
2 

uj k < uj for 7 > — . (45) 
10 



ii) 7 < case 

Through the same calculation, we can show that T_ is positive and T + is negative for all 
negative 7. 



2. 77 > 1 case 

In this case, the behavior of t is singular at r = 0. r and t can be ignored due to the 
divergence of | r From eq.(40) or eq.(41), we know that is always positive definite 
except a singular point r = 0. The condition 77 > 1 is reduced to 

w > - p-iKi-7 2 ) := (46) 

Under this condition, the region of r is divided into —00 < r < and < r < 00. Near 
r = 0, we obtain the behavior of t characterized by the sign of r. When r goes to zero from 
below, t in eq.(41) is written as 

- v (D-iy,B 1 

( K Rh - < 47 » 

When r goes to zero from above, t is reduced as the following 

- Ve (D-lhB 1 

t»~ : r. 48 

(77-I) r^ 1 V ; 

Thus t — > +00 as r — > — but t — > —00 as r — > +0. We also have to examine the behaviors 
of t at r — > ±00. However these have been already described when we discussed the case 
77 < 1. 

To describe the behavior of t as a function of r at r — > ±00 and r — > 0, we classify the 
parameter space of 7 and uj using all results obtained in this section. These are shown in 
Figure 2. 

Figure 2 
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The behavior of the t as a function of r in Figure 2 is summarized as the followings: 

1) In region I, T_ > and T + < (u < ui K ). 

t evolves from finite initial time U to finite final time tf as r runs (— oo, +oo). 

2) In region II, T_ < and T+ < (uj < uj k ). 

t evolves from negative infinity to finite final time tf as r runs (— oo, +oo). 

3) In region III, 71 > and T + < (uj > uj k and uj > uj rj ). 

In this region, because t has a singular behavior at r = 0, the region of r divided into 
— oo < t < and < r < oo. Therefore t has two branches for any given values of 7 and 
uj. For — 00 < r < 0, t evolves from finite initial time U to positive infinity. For < r < 00, 
t evolves from negative infinity to finite final time tf. 

4) In region IV, T_ > and T + > (uj > uo K and uj < u v ). 

In this region, because t has no singularity, t evolves from finite initial time ti to positive 
infinity as r runs (—00, +00). 

5) In region V, T_ > and T + > (uj > uo K and uj > uo rj ). 

By the same reason that explained in region III, t evolves from finite initial time ti to 
positive infinity for — 00 < r < and t evolves from negative infinity to positive infinity for 
< t < 00. 

6) In region VI, T_ < and T+ > (uj < uj k ). 

t evolves from negative infinity to positive infinity as r runs (—00, +00). 

7) In region VII, T_ > and T + > (uj < uj k ). 

t evolves from finite initial time ti to positive infinity as r runs (—00, +00). 

V. THE BEHAVIOR OF THE SCALE FACTOR 

Now we study the behavior of the scale factor a as a function of r. 

A. k < case 

We consider the exponent of the scale factor a(r). Using eq.(10), a(r) is given by 



12 



a(r) 



K 



1 + u;^ — -cr + {1 + u(l - 7)} ln{- cosh(cr)} 
D — 1 c 



5. 



(49) 



In the limit r — > ±00, the scale factor a(r) = e a ^ is rewritten as the following 

a(r) « e H±T , 



(50) 



where _£/± is defined as 



H± = 



2c 



ft 



D-2 r 
l l +u ^- r[ ±{l + u>(l- lf )} 



(51) 



Just as t(r), the behavior of a(r) at r — > ±00 is determined by the sign of H±. Using this 
and the the sign of T±, we can read the behavior of the scale factor a(t) as a function of t 
in the asymptotic regions. 

For negative k {uj < uj k ), H_ > can be written as 



D-2 



(52) 



This means 



1 + ^(1-7) <0 or 

+ >(l+^(l-7)) 5 



(53) 



The first inequality in eq.(53) is equivalent to 



uj < UJ U , 



(54) 



where uj v has been already written in eq.(12). As one can see in Figure 3, uj k < uj u if 7 < jj—^ 
and uj k > uj v if 7 > -^-j-. Together with uj < uj k , the first inequality condition gives the 
region of uj satisfying H > 

1 



uj-i < uj < uj k for 7 < 

The second inequality in eq.(53) is rewritten as 

uj(uj — uj k ) < 0. 
13 



(D-l)- 



(55) 



(56) 



Notice that ui K < if 7 < 2 (d-i) an< ^ > if 7 > 2 (d-i) • Since a; < u; K , we can rewrite 
eq.(56) as the following 

D 



0<u<uj k for 7> 2(D _ 1} . 
As a result, eq.(55) and eq.(57) are the regions of uj satisfying the condition H_ > 0. 



(57) 



Figure 3 



Now we consider the condition H + > 



D-2 

V 1 + ^/J3T > -t 1 + ^( 1 -7)]- 

Like the case if_ > 0, this inequality is divided into two inequalities 



l + w(l-7) >0 or 
2 



(j 1 + LU ^Zjj >(l+w(l-7)) 5 



The first inequality gives the region of u satisfying the condition H + > 

UJ > UJy. 

Using uj > uj-i and k < 0, eq.(60) is rewritten as the following 

1 



c<j_i < uj < uj k for 7 > 



(58) 



(59) 



(60) 



(61) 



D - 1 

The second inequality in eq.(59) has the same region of uj that appeared in eq.(57). Because 
the region of uj in eq.(61) contains the region of uj in eq.(57), eq.(61) is the solution satisfying 
the condition H+ > 0. 



B. k > case 



The exponent of the scale factor a(r) is given by 
2 



K 



l + co^- — \ct -{1 + uj(1- 7)} ln{- I sinh(cr) |} 
JJ — 1 c 



+ B. 



(62) 
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1. r\ < 1 case 



In this case, the scale factor a(r) has no singular behavior. So we investigate the behavior 
of a at r — > ±00. 

In the limit r — > ±00, a(r) is given by 



a t 



(63) 



where if± is defined as 



#4 



2c 



ft 



T {l + U (l-7)} 



The condition ii/l > is exactly equal to eq.(58) due to the sign of ft. When we solve 
eq.(58) under the condition ft > 0, uj > uj k instead of uj < uj k must be applied to the solution. 
Then we obtain the region of uo satisfying the condition > 



uj > uj k for all 7. 

The condition H+ > described by eq.(52) and 00 > uo K gives the region of uo 

D 



(64) 



uj k < uj < for 7 < 



2(D-1)- 



(65) 



2. T] > 1 case 

In this case, we need to investigate the behavior of a(r) at r — > because a(r) has a 
singular behavior at r = 0. In the limit r — > 0, a(r) is written as 

a(r) w e B (g | r |) W . (66) 

For a; > cj^, where —(1 — — uj u ) is negative, a(r) goes to infinite at r — > 0. And for 
c<j < cu^, a(r) goes to zero at r — > 0. The behavior of a(r) at r — > ±00 has been described 
already when we discussed the case 77 < 1. 
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From these studies, we calssify the behavior of a(r) on the parameter space of 7 and uj. 
This is shown in Figure 4. 

Figure 4. 

As shown in Figure 4, we summarize the behavior of a(r) as the followings: 

1) In the region I, f/_ > and H + < (uj < ui K ). 
In the limit r — > ±00, a(r) goes to a zero size. 

2) In the region II, H < and H + > (uj < ui K ). 
In the limit r — > ±00, a(r) goes to an infinite size. 

3) In the region III, H_ > and H + > (uj < ui K ). 

In the limit r — > —00, a(r) goes to a zero size. And in the limit r — > 00, a(r) goes to an 
infinite size. 

4) In the region IV, if_ > and if + > (uj > ui K and uj < uj v ). 

In this region, the behavior of t is not singular. So we need not to consider the behavior of 
a(r) at r = 0. In the limit r — > —00, a(r) goes to a zero size. And in the limit r — > 00, a(r) 
goes to an infinite size. 

5) In the region V, H_ > and H + > (uj > uj k , uj > ui Ti and uj < uj v ). 

In this region, because t has a singular behavior at r = 0, we interpret the behavior of a(r) 
as the following. 

For —00 < r < 0, a(r) goes to a zero size at r — > —00 and r — > 0. For < r < 00, a(r) 
goes to a zero size at r — > and goes to an infinite size at r — > 00. 

6) In the region VI, if_ > and iJ + > (uj > u K , uj > uj v and uj > uj v ). 

For — 00 < r < 0, a(r) goes to a zero size at r — > — 00 and goes to an infinite size at r — > 0. 
For < r < 00, a(r) goes to an infinite size at r — > and at r — > 00. 

7) In the region VII, > and H + < (uj > uj k , uj > ui ri and uj > uj v ). 

For — 00 < t < 0, a(r) goes to a zero size at r — > — 00 and goes to an infinite size at r — > 0. 
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For < r < oo, a(r) goes to an infinite size at r — > and goes to a zero size at r — > oo. 

VI. THE PHASE OF COSMOLOGY 

Using all results obtained from the section IV and V, we now classify the parameter space 
of 7 and oj into several phases and find the behavior of a(t). These phases are characterized 
according to the behavior of a(t). 

Using eq.(19) and eq.(50), in the limit r — > ±oo, a(t) is written as 

a(t) w [T-(t - ti)] H - /T - at -oo, 

a(t) « [T+(t - t f )} H+/T+ at t -> oo, (67) 

where tj and t/, which are defined in section IV, are real constant. Notice that ti (tf) 
becomes the starting point (the ending point) in the case T_ > (T+ < 0) and that ti (tf) 
can be neglected in the case T_ < (T + > 0) because t — > ±oo as r — > ±oo. 
Now, we explain two examples: 

i) For T_ < and iJ_/T_ > 0, T_(t — ti) is positive and a(t) goes to positive infinite at 
t — > — oo. 

ii) For T_ > and iJ_/T_ > 0, i can be defined in the region t > U only. So £ — tj is 
positive and a(t) goes to zero at t — > t^. 

Other cases can be analysed through the same method. 

In the region u > u K and r] > 1, we must investigate the behavior of a(t) at r — > 0. From 
eq.(47), eq.(48) and eq.(66), a(t) is obtained as the following 

a(t) ss -Ex | t | S=TTra , (68) 

where 

2(l- 7 )(^-^) R r. 2(0-1)7(1-7)^-^) 1 

E = [g(^ - l)] (v-i)M e L PFTira J 
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is a positive value because q and (1 — 7) are positive in the previous definition. a(t) goes to 
zero at t — > ±00 (r — > ±0) in the case a; < u;„ and a(t) goes to infinite at t — > ±00 in the 
case uj > uj v . 

Figure 5 

As shown in Figure 5, using the sign of T± and H± with the consideration of the behavior 
of a(t) at t — > 0, the behavior of a(t) in each region is characterized as the followings: 

1) In region I, T_ > 0, T+ < 0, H > 0, and H + < 0. 

The universe evolves from a zero size at finite initial time ti to a zero size at finite final time 

tf- 

2) In region II, T_ < 0, T + < 0, if_ > 0, and H+ < 0. 

The universe evolves from a zero size at negative infinity to a zero size at finite final time 

tf- 

3) In region III, T_ > 0, T + > 0, H_ > 0, and H + > 0. 

The universe evolves from a zero size at finite initial time ti to an infinite size at positive 
infinity. 

4) In region IV, T_ < 0, T + > 0, if_ < 0, and #+ > 0. 

The universe evolves from an infinite size at negative infinity to an infinite size at positive 
infinity. 

5) In region V, T_ > 0, T + > 0, H < 0, and H + > 0. 

The universe evolves from an infinite size at finite initial time ti to an infinite size at positive 
infinity. 

6) In region VI, T_ > 0, T + > 0, H_ > 0, and H + > 0. 

The universe evolves from a zero size at finite initial time ti to an infinite size at positive 
infinity. 

From 7) to 11), we consider the case 77 > 1 and u > uo K in which t has a singular behavior 
at r — > 0. In these cases, we can divide the region of r into — 00 < r < and < r < 00. 

18 



Therefore we obtain two branches of a(t) having the different behaviors in each region of r. 

7) In region VII, T_ > 0, T + > 0, > 0, and H + > 0. 

In the region — oo < r < 0, the universe evolves from a zero size at finite initial time U 
to a zero size at positive infinity. 

In the region < r < oo, the universe evolves from a zero size at negative infinity to an 
infinite size at positive infinity. 

8) In region VIII, T_ > 0, T + < 0, > 0, and H + > 0. 

In the region — oo < r < 0, the universe evolves from a zero size at finite initial time U 
to a zero size at positive infinity. 

In the region < r < oo, the universe evolves from a zero size at negative infinity to an 
infinite size at finite final time tf. 

9) In region IX, T_ > 0, T + > 0, H > 0, and H + > 0. 

In the region — oo < r < 0, the universe evolves from a zero size at finite initial time U 
to an infinite size at positive infinity. 

In the region < r < oo, the universe evolves from an infinite size at negative infinity 
to an infinite size at positive infinity. 

10) In region X, T_ > 0, T + < 0, H > 0, and H + > 0. 

In the region — oo < r < 0, the universe evolves from a zero size at finite initial time ti 
to an infinite size at positive infinity. 

In the region < r < oo, the universe evolves from an infinite size at negative infinity 
to an infinite size at finite final time tf. 

11) In region XI, T_ > 0, T + < 0, H > 0, and H + < 0. 

In the region — oo < r < 0, the universe evolves from a zero size at finite initial time U 
to an infinite size at positive infinity. 

In the region < r < oo, the universe evolves from an infinite size at negative infinity 
to a zero size at finite final time tf. 
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VII. DISCUSSION AND CONCLUSION 



In this paper we studied the effect of the gas of solitonic p-brane by treating them as 
a perfect fluid in the Brans-Dicke theory. We found exact cosmological solutions for any 
Brans-Dicke parameter uj and for general constant 7 and classified the cosmology of the 
solutions according to the parameters involved. We assumed that the universe is dominated 
by one kind of p-brane and they are treated as perfect fluid. We found the analytic solution 
which is singularity free for some 7 and uj. It is very interesting that a(t) has no initial 
and final singularities at finite initial and final cosmic time in region IV and in region VII. 
a(t) has also an inflation behavior in region VII. So we need to study more intensively the 
behavior of a(t) in these regions. 

Presumably the value of 7 as well as uj should be fixed once p is fixed. Without knowing 
the value of 7 for a given p, the classification was the best thing we could do. It would be 
very interesting to determine the parameter 7 for the given p. Also we need more rigorous 
justification of our basis for the p-brane cosmology. If what we took as basis goes wrong, 
then what we have done is just Brans-Dicke cosmology in the presence of some perfect fluid 
type matter. We wish that more study of the effect of the solitons in the string cosmology 
be done in the future. 
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Caption 



Figure 1. In 4 dimension (D = 4), all functions (u K ,u v , ■ ■ ■), defined by the relation 
between t and r, are presented on a parameter space of 7 and ui. 

Figure 2. In 4 dimension (D = 4), the parameter space is classified by the relation 
between t and r. 

Figure 3. In 4 dimension (D = 4), all functions (cu K , cu u , ■ ■ •), defined by the relation 
between a(r) and r, are presented on a parameter space of 7 and cj. 

Figure 4. In 4 dimension (D = 4), the parameter space is classified by the behavior of 
a(r). 

Figure 5. In 4 dimension (D — 4), the parameter space is classified by the behavior of 
a(t). 
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